Entanglement dynamics of a strongly driven trapped atom 
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We study the entanglement between the internal electronic and the external vibrational degrees 
of freedom of a trapped atom which is driven by two lasers into electromagnetically-induced trans- 
parency. It is shown that basic features of the intricate entanglement dynamics can be traced to 
Landau-Zener splittings (avoided crossings) in the spectrum of the atom-laser field Hamiltonian. 
We further construct an effective Hamiltonian that describes the behavior of entanglement under 
dissipation induced by spontaneous emission processes. The proposed approach is applicable to a 
broad range of scenarios for the control of entanglement between electronic and translational degrees 
of freedom of trapped atoms through suitable laser fields. 
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A trapped atom driven by two laser beams under 
conditions of electromagnetically-induced transparency 
(EIT) represents a realistic model for current experi- 
ments on ion or atom trapping [TJ [5] . Quantum correla- 
tions between the translational and the electronic degree 
of freedom of the trapped atom play a crucial role in the 
cooling dynamics of trapped atoms and ions. Indeed, 
such correlations are responsible for the transfer of vi- 
brational energy of the trapped atom into excess energy 
of the scattered radiation field [3]. 

Here, we explore the dynamical behavior of entangle- 
ment [H [5] between the external vibrational degree of 
freedom of the trapped atom and its internal electronic 
degree of freedom during the cooling process. Solving 
numerically the master equation [6 which describes the 
dynamics of the composite state p(t) of the trapped atom- 
laser field system, we investigate the evolution of quan- 
tum correlations by use of the negativity [7] 

N{t)=N{p(t))= l -(\\p T {t)\\-l), (1) 

where p T (t) represents the partial transpose of pit), and 
|| • || denotes the trace norm. The negativity M{p) is a 
nonnegative function which quantifies the degree of en- 
tanglement in a mixed quantum state p. It vanishes for 
separable, classically correlated states, and takes on a 
maximum value if p is a pure, maximally entangled state. 

In the trapped atom system the Lamb-Dicke parameter 
r\ = fc/v / 2mw (H=l) controls the probability of changing 
the vibrational quantum number n in electronic transi- 
tions (5]. For small r\ the photon-recoil energy k 2 /2m 
is an inconsiderable fraction of the trap vibrational fre- 
quency lo. In the Lamb-Dicke regime, r\\Jn + 1 <C 1, only 
the first sideband transitions, An = ±1, are significant, 
their probability being proportional to (n + l)?7 2 , where 
n is the vibrational quantum number. Even outside the 
Lamb-Dicke regime a rapid removal of vibrational energy 
of a trapped atom occurs under suitable experimental 
parameters, i.e., a difference of ac-Stark shifts of the two 
ground states of the order of w [3j [6j [9] . Under these con- 
ditions, with blue detuning of two counter-propagating 



lasers, a vibrationally excited atom can be driven into 
states with mean values (n) near zero. 

An intricate behavior of the negativity is observed dur- 
ing such a vibrational cooling process, as is shown in 
Fig. [I] We see a rapid build-up of entanglement, a sub- 
sequent exponentially damped oscillation, followed by a 
precipitous drop and a power law behavior and, finally, 
a rebirth of entanglement towards the non-equilibrium 
stationary state. 
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FIG. 1: Dynamics of the negativity during the cooling process 
for the initial vibrational level n< = 10 for three values of the 
excited state width Y and r\ = 0.1. The top scale gives the 
mean vibrational quantum number for the 6 MHz result. The 
dotted line represents a power law decay, M oc £~ 3//2 . 

We demonstrate that this complex behavior of the en- 
tanglement dynamics can be related to Landau-Zener 
splittings or avoided crossings in the spectrum of the 
atom-laser field Hamiltonian. The dominant source of en- 
tanglement is thus a near degeneracy of vibronic states, 
which appears as a consequence of the ac-Stark shift, 
suitably chosen for EIT-cooling. Two-photon transitions 
associated with a change in vibrational quantum number 
occur and lead to a periodic oscillation of entanglement. 
In addition, we derive an effective Hamiltonian that al- 



2 



lows the modeling of the negativity decay in the presence 
of dissipation induced by spontaneous emission processes. 
Our model is based on the quantum master equation 



dt 



p(t) = -i[H\p(t)]+£p{t) 



(2) 



for the interaction picture density matrix pit) represent- 
ing the combined state of vibrational and electronic de- 
grees of freedom of the trapped atom. The total Hamil- 
tonian 



(3) 



of the model consists of three parts: H cm = Loa^a de- 
scribes the vibrational degree of freedom of a harmoni- 
cally trapped atom with raising and lowering operators 
of and a, and H c \ = A(|l)(l| + |2)(2|) represents the elec- 
tronic degree of freedom with detuning A. The electronic 
states are denoted by i = 1,2,3, and form a A-type 
level structure sketched in Fig. [2^,. The Hamiltonian 
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FIG. 2: a) Trapped 3-level atom, b) energies of dressed vi- 
bronic states as a function of the normalized ac-Stark shift. 
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|3)(2|+h.c. (4) 



models the interaction with two counter-propagating 
laser beams, the most efficient configuration for cooling, 
where x denotes the atomic position and k the wavenum- 
ber (ki m k 2 = k). Note that H^ nt depends on the Lamb- 
Dicke parameter 77 through kx = rj{a) + a). Finally, dis- 
sipation caused by spontaneous emissions is modelled by 
the Lindblad superoperator (see, e.g., Ref. [TU] ) 



£ '=EE 



j=l,2q=± 



(5) 



where T denotes the spontaneous emission rate and 
a j+ = \j){i\e lkx , (j j— = \j){2,\er lkx . The Lindblad oper- 
ators <Tj_|_ describe the transitions from the excited elec- 
tronic state 1 3) to the ground states \j = 1, 2) with a pho- 
ton emitted parallel to the harmonic oscillator axis, while 



the <7j- describe the transitions from |3) to \j — 1,2) 
with a photon emitted counter-propagating to the har- 
monic oscillator axis. This is a good compromise short 
of doing a full 3D simulation as the recoil of the transla- 
tions! motion is zero for emissions perpendicular to the 
trap axis, while it leads to Dopplcr heating or cooling for 
emissions along the trap axis. Typical results obtained 
by numerically solving Eq. ^ [3l [6] are given in Fig. [I] 
We start our discussion of the entanglement dynam- 
ics by considering first the case T = 0, neglecting com- 
pletely effects from spontaneous emission processes. The 
dynamics is then given by the Schrodinger equation 



dt 



m)) = -^\m) 



(6) 



Our goal is to describe the basic features of the dynamical 
behavior of the negativity in terms of the properties of 
the spectrum of the full Hamiltonian H v given by Eq. Q . 
For vanishing Lamb-Dicke parameter, rj — 0, the external 
and the internal degrees of freedom of the atom decouple, 
such that the eigenstates of H° = H cm + H c \ + H? nt are 
given by tensor product states \ipi,n) = \ipi) (g> \n) with 
corresponding energy eigenvalues, e^ n = £, + ncu. Here, 
n = 0, 1, 2, . . . denotes the vibrational quantum number. 
The \ifii) are dressed electronic states which are defined as 
normalized eigenstates of H c \ + H® nt with the eigenvalues 



A- 



£3 



sac- 



(7) 



9l 



9l 



A ] is the ac-Stark 



where e A c = \ (V^ 1 
shift. Note that we have labelled the dressed states such 
that for gi 3> <?2 the state \ifi) has the main weight on 
the electronic state \i), and that ^2) represents the dark 
state which decouples from the laser fields for r; = 0. 
It follows from the above equations that the levels £2, ra 
and £i in _i cross each other (see Fig. [2Jd) if the resonance 
condition 



£ac 



4A 



(8) 



is satisfied. The approximation in Eq. ^ refers to the 
case gi t 2 <§C A which is assumed throughout the paper. 

For non-vanishing Lamb-Dicke parameter, 77 > 0, the 
electronic and the translational degrees of freedom of the 
atom are coupled. This coupling lifts the degeneracy of 
the levels £2,n and £i in _i at the crossing point to yield an 
avoided crossing or Landau-Zener splitting [TTHT3] with 
a certain energy gap AE n , as is illustrated in Fig. [2}d. 

Employing degenerate perturbation theory it can be 
shown that the eigenstates at the center of the avoided 
crossing are given by the even and odd linear combina- 
tions of the unperturbed (crossing) states. More pre- 
cisely, this means that with an appropriate choice of the 
phase of the states \<p2,n) and \ipi,n— 1) the eigenstates 
of the full Hamiltonian at resonance have the approxi- 
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mate form 



(\ 



\\ )=^=\^ipi)9\n)±\ Vl )2)\n-l)) (9) 



with the corresponding energy eigenvalues 
E± = A + nu± -. 



(10) 



Considering as initial state |^(0)) = \cp2,n), it follows 
that the solution of the Schrodinger equation is given by 



E+t\ 



1X4 



e-^-*|X-! 



-iE+t 



1 



( e -^+< - < 



The Schmidt coefficients 
1 , 



e- i£ - 4 ) fa) ® \n) 

' E - t )\ Vl )®\n-l). (11) 
of this state take the form 



"1.2 



which immediately leads to the negativity: 



1 



M{t)=a x a 2 = ~\sm(AE n -t)\ 



(12) 



(13) 



The same formula is obtained if one considers the initial 
state |^(0)) = \(fi,n — 1). Thus, we see that the en- 
tanglement dynamics is determined by a single spectral 
quantity, namely the energy gap AE n , describing an os- 
cillation of the negativity with the period T = ir/AE n . 
In Fig. [3] we compare the negativity obtained from the 
numerical solution of the full Schrodinger equation §6§ 
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FIG. 3: Time dependence of the negativity obtained from 
the Schrodinger equation |6| for initial vibrational levels 
rii = 3 and 30. Parameters: 77=0.1, a>=0.03, A=15, <7i = 1.34, 
ffi/<?2=10 (all frequencies are given in units of 27rMHz). The 
broken lines represent the prediction according to Eq. (131. 



with the formula ( 13 ) for various initial states, which 



clearly demonstrates the accuracy of the approximation 
based on a single Landau-Zener state pair ®. The nega- 
tivity predicted by Eq. ^ slightly exceeds the value 1 /2 
as it also accounts for the small contributions to entan- 
glement by the excited state levels \<p3,n). We empha- 



size that the argument which leads to the result ( 13 1 is 



quite general and should be applicable to many differ- 
ent physical scenarios: It only requires the existence of 
a Landau-Zener splitting between tensor product states 
on some composite Hilbert space. Moreover, depending 



only on the energy gap, formula ( 13 ) can be applied even 



if the structure of the factors of the tensor products is 
not known. 

Next we discuss the influence of dissipation on the en- 
tanglement dynamics which is described by the Lindblad 
dissipator of Eq. ([5| . Figure [I] shows examples for the 
negativity dynamics obtained from numerical solutions 
of the full master equation ([2| . We see that the negativ- 
ity follows a damped oscillation for small T, where the 
revivals of the negativity come in pairs of nearly equal 
heights. At a certain threshold value of T the oscillations 
die out and a monotonic decay of the negativity after the 
first bump is observed. 

The entanglement behavior under dissipation can be 
modelled through the Schrodinger equation 



dt 



\m) = 



with the effective Hamiltonian 



i?|3)<3| 



(14) 



(15) 



This is a non-Hermitian Hamiltonian derived from the 
quantum master equation ([2]): It describes the coherent 
motion generated by H v as well as the loss term (anti- 
commutator) of the dissipator |5]) which leads to finite 
widths of the dressed states of the atom-field system 
given by the imaginary parts of the complex eigenval- 
ues. Employing the effective Schrodinger equation (14) 
we thus neglect the gain terms of the dissipator ([5]). 
These latter terms describe real emission processes and 
induce, at least partially, transitions into other Landau- 
Zener state pairs which add incoherently and, thus, lead 
only to negligible contributions to the negativity. 

In the framework of this approximation we can extend 
the Landau-Zener state pair calculation to the effective 
Hamiltonian (15). Under the conditions gi 2 <C A and 



I< A the eigenvalues of H^ s are 



£■2 



£3 



= A 
= A, 



9l 



92 



■ 7i 



9l 



4A 



9l 



4A 
-.9l)/(4A 2 ). 



(16) 
(17) 

7i), (18) 
We denote the eigenstates 



2 V 



where 71 = ■ 
corresponding to the eigenvalues £1,2 again by |y>i,2) 5 an d 
the Landau-Zener pair by \(p\,n— 1) and \tf2,n). Using 



these states as basis states, the effective Hamiltonian ( 15 ) 



reduces within the two-state resonance approximation to 
the matrix 



A + nw-i7i/2 AE n /2 
AE n /2 A + nu 



(19) 
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This Hamiltonian may be viewed as describing damped 
Rabi-type oscillations between the Landau-Zcner pair 
states with Rabi frequency AE n , where only the state 
\ifi,n — 1) is damped, since the width of the dark state 
\(f2,n) is practically zero. With the effective Hamilto- 
nian dl9|) it is now easy to determine the solution of the 



Schrodinger equation ( 14 1 corresponding to the initial 
state IV'(O)) = \ip2,n), and to derive for the negativity, 



Af(t) 



AE n 



sin vt 



_7i 
2v 



2 (1 — COSVt) 



- 7l t/2 



(20) 



jf /4. As is illustrated in Fig. [4] 



where v = \J (AE n ) 2 

the simple formula |20j provides an excellent approxi- 
mation of the entanglement dynamics for 71/2 < AE n , 
and even for intermediate values of the spontaneous emis- 
sion rate the qualitative behavior of the negativity is 
reproduced. It is remarkable that Eq. (201 correctly 
describes the transition from the underdamped to the 
overdamped motion of the negativity at 71/2 = AE n , 
compare Figs. [3J; and [4ji. Using the approximations 
AE n s» rjy/ngi g 2 / A and 71 ~Tg 2 /{AA 2 ) this transition is 
predicted to occur at the value T = 8ri\/nAg2/gi, which 
nicely fits to the results of our numerics. 
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FIG. 4: The negativity obtained from the quantum master 
equation Q for various F, as also shown in Fig. [l] Fits of 
Eq. ( |20[) are shown by the broken lines. Conditions as given 
in Fiiq| 

The long-time behavior of the negativity for larger val- 
ues of the spontaneous emission rate T exhibits an in- 
teresting and complex structure, see Fig. [T] The expo- 
nential decay of the initial peak of the negativity is fol- 
lowed by a precipitous drop until the rate of decay slows 
down and a smooth transition to an algebraic decay of 
the form Af(t) cx t~ z / 2 emerges, see the broken line in 
Fig. [T] While a detailed discussion is beyond the scope 
of this note, we mention that the Landau-Zener picture 
can also offer an explanation here. A power law gov- 



erns the negativity dynamics until a rebirth of entangle- 
ment towards the final value N{oo) = A/"(p s tat) occurs. 
The final value of the negativity represents the entan- 
glement in the unique non-equilibrium stationary state 
Pstat of the master equation During the cooling the 
atom is driven into a mixed stationary state which has a 
dominant contribution from the state \ip2,0), the prod- 
uct of the dark state and the vibrational ground state 
of the trap. However, p s tat contains a small amount of 
entanglement which is controlled by the small steady- 
state heating caused by off-diagonal transitions from the 
dark state. Also here the LZ-model offers a quantitative 
explanation. The interaction of the lowest dark state, 
(</?2, Ol-ffjnJ^i, 1) = 77^1(72/ (4 A), leads to a steady state 
negativity, Af (oo) =77 .92/91 for T— >0, in agreement with 
the predictions from Eq. |2]) for the examples in Fig. [I] 

The Landau-Zener mechanism suggest various meth- 
ods for the control of entanglement between electronic 
and translational degrees of freedom of trapped atoms 
by means of suitable laser fields. One possibility is the 
switching of fields at appropriate points of time, e.g., at 
the instant of maximal entanglement. Since the mani- 
fold spanned by the Landau-Zener state pair has only a 
negligible contribution from the excited electronic state, 
the negativity practically freezes upon rapidly switching- 
off the lasers. This method can thus create permanent 
entanglement of the trapped atom. A further possibil- 
ity is to use suitable laser pulses in order to generate 
entanglement through the control of adiabatic/diabatic 
Landau-Zener transitions. A detailed discussion of these 
points will be given in a forthcoming paper. 
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